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^ ■ Abstract 

Q\^ [ The generalized sine-Gordon (sG) equation utx = (1 + ^d^) sin u was derived as an inte- 

grable generalization of the sG equation. In a previous paper (Matsuno Y 2010 J. Phys. 

^ I A: Math. Theor. 43 105204) which is referred to as I, we developed a systematic method 

for solving the generalized sG equation with z/ = — 1. Here, we address the equation with 
u = 1. By solving the equation analytically, we find that the structure of solutions differs 
substantially from that of the former equation. In particular, we show that the equation 

VD ' exhibits kink and breather solutions and does not admit multi-valued solutions like loop 

lO ■ solitons as obtained in I. We also demonstrate that the equation reduces to the short pulse 



and sG equations in appropriate scaling limits. The limiting forms of the multisoliton so- 
lutions are also presented. Last, we provide a recipe for deriving an infinite number of 
conservation laws by using a novel Backlund transformation connecting solutions of the 
sG and generalized sG equations. 
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1. Introduction 

The generalized sine-Gordon (sG) equation 

Utx = (1 + ^d"^) sin u, (1.1) 

where u = u{x,t) is a scalar-valued function, z/ is a real parameter, d^ = d'^/dx^ and the 
subscripts t and x appended to u denote partial differentiation, has been derived in [1] 
using bi-Hamiltonian methods. In the case of z/ = —1, its integrability was established by 
constructing a Lax pair associated with it and the initial value problem was formulated 
for decaying initial data by means of the inverse scattering method [2]. Quite recently, 
we developed a systematic method for solving equation (1.1) with u = —1 and obtained 
soliton solutions in the form of parametric representation [3] , which will be referred to as 
I hereafter. We showed that the solutions exhibit various new features which have never 
seen in the sG solitons. We also demonstrated that the generalized sG equation reduces 
to the physically important short pulse and sG equations in appropriate scaling limits. 
In this paper, we consider equation (1.1) with u = 1 

ut, = {l + dD sin u. (1.2) 

Despite its resemblance to equation (1.1) with u = —1, the analysis of equation (1.2) by the 
inverse scattering method has not been done and explicit solutions are still unavailable [2] . 
This is a motivation why we address equation (1.2). We use an exact method of solution 
developed in I to construct soliton solutions and investigate their properties. One of the 
remarkable features of equation (1.2) is that it does not admit multi- valued solutions like 
loop solitons as obtained in I. It is therefore of considerable interest to study the structure 
of solutions in conparison with those of equation (1.1) with u = —1. 

This paper is organized as follows. In section 2, we summarize a direct method of 
solution in such a manner that is relevant to equation (1.2). Specifically, a hodograph 
transformation makes it possible to transform equation (1.2) into an integrable system of 
nonlinear partial differential equations (PDEs). This system is further recast to a system 
of bilinear equations through the dependent variable transformations. The solutions to 
the bilinear equations can be constructed following the standard procedure in the bilinear 



transformation method [4, 5]. Inverting the relationship that determines the hodograph 
transformation, we obtain a parametric representation of the solution of equation (1.2) 
in terms of tau functions. In section 3, we present soliton solutions and investigate 
their properties focusing mainly on 1- and 2-soliton and general multisoliton solutions. 
Throughout this paper, we use the term 'soliton' as a generic name of elementary solutions 
such as kink, antikink and breather solutions. In section 4, we show that equation (1.2) 
reduces to the short pulse equation in an appropriate scaling limit. The limiting form 
of the multisoliton solution is also presented. We find that the short pulse equation 
exhibits a novel type of singular soliton solutions. Last, we briefly discuss the reduction 
to the sG equation, reproducing the known results about the multisoliton solution and 
its characteristics. In section 5, we provide a recipe for deriving an inflnite number of 
conservation laws starting from those of the sG equation. It is based on a novel Backlund 
transformation connecting solutions of the sG and generalized sG equations. Section 
6 is devoted to conclusion where we discuss some open problems associated with the 
generalized sG equation. In Appendix A, we present a method for obtaining solutions of 
traveling type and recover the 1-soliton solutions. 

2. Exact method of solution 

2.1. Hodograph transformation 

Here, we summarize an exact method of solution which consists of a sequence of nonlinear 
transformations. First, we introduce the new dependent variable r in accordance with the 
relation 

r'^ = l-ul, (0<r<l), (2.1) 

to transform equation (1.2) into the conservation law of the form 

rt — (r cos u)x = 0. (2.2) 

This expression makes it possible to deflne the hodograph transformation {x, t) — )■ {y, r) 
by 

dy = rdx + r cos u dt, dr = dt. (2.3) 



The X and t derivatives are then rewritten in terms of the y and r derivatives as 

d d d d d ,^ ^, 

ox ay at or oy 

With the new variables y and r, (2.1) and (2.2) are recast into the form 

r2 = 1 - r'^ul, (2.5) 

fij +(cosm)j, = 0, (2.6) 

respectively. Further reduction is possible if one defines the variable by 

Uy = sinh (f), = (f){y, r). (2.7) 

It follows from (2.5) and (2.7) that 

- = cosh 0. (2.8) 

r 

Substituting (2.7) and (2.8) into equation (2.6), we find 

(pr = sin u. (2.9) 

If we eliminate the variable from (2.7) and (2.9), we obtain a single PDE for u 

"""^ - sin u. (2.10) 



nr+l?y 



Similarly, elimination of the variable u gives a single PDE for 



^ry 



sinh 0. (2.11) 



By inverting the hodograph transformation (2.4) and using (2.8), the equation that 
determines the inverse mapping {y, r) — )> {x, t) is found to be governed by the system of 
linear PDEs for x = x{y, r) 

Xy = cosh 0, (2.12a) 

Xt = — cos u. (2.126) 

It is important that the integrability of the system of equations (2.12) is assured by (2.7) 
and (2.9). Thus, given and u, the above system of equations can be integrated to 



give the parametric expression of x in terms of y and r. Remarkably, we were able to 
perform the integration analytically. This last step for constructing solutions is the core 
of the present analysis. The main results associated with soliton solutions will be given 
by theorem 2.1 and theorem 2.2 below. 

2.2. Bilinear formalism 

Here, we develop a method for solving a system of PDEs (2.7) and (2.9). Let a and a' be 
solutions of the sG equation 

a-ry = sin a, a = a{y, r), (2.13a) 

aly = sma', a' = a'{y,T). (2.136) 

The solutions of the above equations can be put into the form [6-9] 

f 
(T = 2ilny, (2.14a) 

q' 
cr' = 2iln-. (2.146) 

9 
For soliton solutions, the tau functions /, /', g and g' satisfy the following system of 

bilinear equations [6]: 

D,DJ-f = ^{f-f"), (2.15a) 

DrDyf'-f' = ^if"-f), (2.156) 

D,Dyg-g = ^{g^-g''), (2.16a) 

D,Dyg'-g' = ^{g"-g'), (2.166) 

where the bilinear operators D^- and Dy are defined by 

D^D^f ■g = {dr- dr'T (dy - dyrf{r,y)9{r',y%'=r,y'=y, {m,n = 0, 1,2, ...). 

(2.17) 
Now, we seek solutions of equations (2.7) and (2.9) of the form 

F' 
M = iln— , (2.18a) 

F 



= ln^, (2.186) 

where F, F\ G and G' are new tau functions. If we impose the condition 

F'F = G'G, (2.19) 

among these tau functions, then equations (2.7) and (2.9) can be transformed to the 
following bilinear equations 

iDyF' ■F=-{G'^-G^), (2.20) 

iD,G'-G=-{F^-F'^), (2.21) 

respectively. The proposition below provides the tau functions F, F', G and G" in terms 

of /, /', 9 and g'. 

Proposition 2.1. // we impose the conditions for the tau functions f, f',g and g' 

^Dyf-g'=^ifg'-f'g), (2.22a) 

iDyf'-g = ^{f'g-fg'), (2.226) 

iDJ-g = -^ifg-f'g'), (2.23a) 

iD,f-g' = -^{fg'-fg), (2.236) 

then the solutions of bilinear equations (2.20) and (2.21) subjected to the condition (2.19) 
are given by 

F = fg, F' = f'g', (2.24a) 

G = fg', G' = f'g. (2.246) 

Proof. It is obvious that the tau functions (2.24) satisfy the condition (2.19). We first 

prove (2.20). Substituting (2.24) into the left-hand side of (2.20) and using (2.22), we 
find that 

iDyF' ■ F = i{{DJ' ■ g)fg' - {DJ ■ g')fg} 
= \{G"-G% 



which is (2.20). The proof of (2.21) can be done in the same way by using (2.23). D 

2.3. Parametric representation 

We demonstrate that the solution of equation (1.2) admits a parametric representation. 
The following relation is crucial to integrate (2.12): 

Proposition 2.2. cosh is given in terms of the tau functions /, f ,g and g' as 

cosh 0=1 + i (in^] . (2.25) 

Proof. Using (2.22), one obtains 



,nn'^\ ^_,MljL.,D.r-g 



f'fj, f'J f's 

1 (fsT + (far _ ^ p^26) 



(2.27) 



2 ffg'g 

On the other hand, it follows from (2.18b) and (2.24b) that 

cosh0 = -^- + - 

_ 1 {fg'f + {fg? 
2 ffg'g 

The relation (2.25) follows immediately by comparing (2.26) and (2.27). D 

Integrating (2.12a) with (2.25) by y yields the expression of x 

x = y + i\nj^ + d{T), (2.28) 

where d is an integration constant which depends generally on r. The expression (2.28) 
now leads to our main result: 

Theorem 2.1. The solution of equation (1.2) can he expressed by the parametric repre- 
sentation 

u{y,T) = i\n^, (2.29a) 

x{y, T)=y-r + ilnj^ + yo, (2.296) 



where the tau functions f,f',g and g' satisfy equations (2.15), (2.16), (2.22) and (2.23) 
and yo is an arbitrary constant independent of y and r. 

Proof. The expression (2.29a) for -u is a consequence of (2.18a) and (2.24a). To prove 
(2.29b), we substitute (2.28) into (2.12b) and obtain the relation 

i Mn|;^j + d'{T) = - cos u. (2.30) 

The left-hand side of (2.30) is modified by using (2.23) whereas the right-hand side can 
be expressed by /, f',g and g' in view of (2.29a). After a few calculations, we find that 
most terms are cancelled, leaving the equation d'{T) = —1. Integrating this equation, one 
obtains (i(r) = —r + yo, which, substituted into (2.28), gives the expression (2.29b) for x. 

D 

An interesting feature of the parametric solution (2.29) is that it never exhibits singu- 
larities as encountered in the case of equation (1.1) with u = —1 (see I). To demonstrate 
this, we calculate u^ from (2.7) and (2.8) and obtain 

Ux = ruy = tanh 0. (2-31) 

The relation (2.31) implies that u^ always takes finite value and singular solutions such 
as loop solitons and multi- valued kinks obtained in I never exist for equation (1.2). See 
Figures 2-4 in I. It should be remarked, however, that the above relation does not exclude 
the presence of singular solutions like peaked waves (or peakons), for instance which have 
a finite discontinuity in their slope at the crest. The construction of the latter type of 
solutions is of some interest from the mathematical point of view. 

2.4. Multisoliton solutions 

The last step in constructing solutions is to find the tau-functions /, f',g and g' for the 
sG equation which satisfy simultaneously the bilinear equations (2.22) and (2.23). The 
following theorem establishes this purpose: 

Theorem 2.2. The tau-functions f, f',g and g' given below satisfy both the bilinear forms 



(2.15) and (2.16) of the gG equation and the bilinear equations (2.22) and (2.23), 



f = J^ exp 

11=0,1 



■ N 

J2 f^j (O + c^j + ^ i) + J2 /^i^fc^jA 

.j=l l<j<k<N 



/' = IJ exp 

M=0,1 


' N 

J2 N (ij + c?, - 1 i) + Y^ fijfikijk 

.j=l l<j<k<N 


c, = ^ exp 

M=0,1 


' N 

.3=1 i<j<k<N 


g' = "^ exp 

fi=0,l 


' N 

J2 f^j (o - ^i - f i) + Yl ^^3^^kljk 

_j=l l<j<k<N 



where 



1 

4j = PjV +—r + 4iO, 

Pj 



oHik 



Pj - Pk 
Pj + Pk 



ipj 



(j = l,2,...,iV), 
{j,k = l,2,...,N-jy^k), 
(j = l,2,...,iV). 



(2.32a) 

(2.326) 
(2.33a) 
(2.336) 

(2.34a) 
(2.346) 

(2.34c) 



1 - ipj 

Here, pj and $^jq are arbitrary complex parameters satisfying the conditions pj ^ ±pk for 
j ^ k, i = \/—l and N is an arbitrary positive integer. The notation X]u=o i 'i'^^pli^s the 
summation over all possible combination of fii = 0, l,/i2 = 0, 1, ...,/iAf = 0, 1. 

Proof. It has been shown that / and /' from (2.32) and g and g' from (2.33) satisfy 
the bihnear equations (2.15) and (2.16), respectively [6]. To prove that (2.32) and (2.33) 
satisfy the bihnear equations (2.22) and (2.23), we replace the parameters pj by —ipj {j = 
1, 2, ..., A^) and shift the phase parameters as ^jo — )■ ^jo — dj {j = 1,2,..., A^) and at the 
same time change the variables y and r as y — )■ iy, r — )■ — ir, respectively. It then turns 
out that ^j and 7^^ defind respectively by (2.34a) and (2.34b) remain the same form. 
Thus, the proof reduces to the theorem 2.2 in I. See also Remark 2.5 in I. D 

The parametric solution (2.29) with (2.32) and (2.33) is characterized by the 2N 
complex parameters pj and C,jo {j = 1,2,..., A^). It produces in general the complex- valued 
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solutions. The real-valued solutions are obtainable if one imposes certain conditions on 
these parameters. Actually, there arise various type of solutions depending on values of 
the parameters. These solutions include kinks, antikinks and breathers. Among them, we 
consider following three types: 

Type 1: Kink solution 

First, let pj and ^jq {j = 1,2, ..., A^) be real quantities. It then follows from (2.34c) that 
d* = —dj {j = 1,2,. ..,N) where the asterisk denotes the complex conjugate. Taking 
account this fact in (2.32) and (2.33), we find that /' = g* and g' = /*. These relations 
are substituted into (2.29) to give the real parametric solution 

f*g* 

u{y,T) = i\n— — , (2.35a) 

xiy,T) =y-T + i\n— + yo. (2.356) 

As will be demonstrated in section 3, this solution yields the A^-kink solution if the 
parameters pj are all positive. The expression u^ represents the A^-soliton solution which 
decays at infinity. Remarkably, the amplitude of each soliton does not exceed 1 by virtue of 
the relation (2.31). The A^-antikink solution is also available if these parameters are set to 
be negative. We can also produce kink-antikink solutions. For example, let A^ = M + M' 
where M and M' are positive integers and assign M positive and M' negative values for 
A^ parameters pj {j = 1,2, ...,N). Then, the parametric solution represents the solution 
describing the interaction among M kinks and M' antikinks. 

Type 2: Breather solution 

The second example is the parameterization which gives rise to the breather solution. To 
show this, we put A^ = 2M where M is a positive integer, and specify the parameters pj 
and^j-o (j = 1,2, ...,2M) as 

P2J-1 = P*2j, 6,-1,0 = C2j,o, (j = 1, 2, ..., M). (2.36) 

We see from (2.34) and (2.36) that 6i-i = ^2j^ c?2j-i = -d^j (j = 1,2, ...,M). It turns 
out from (2.32) and (2.33) that /' = g* and g' = /*. Then, the solution can be written 
in the same form as (2.35). 

10 



Type 3: Kink-breather solution 

Let A^ = 2M + M' where M and M' are positive integers. In addition to the parameteriza- 
tion given by (2.36), the 2M' parameters Pj(> 0) and ^^o (j = 2M+1, 2M+2, ..., 2M+M') 
are chosen to be reaL Then, the parameteric solution (2.35) represents the solution de- 
scribing the interaction among M breathers and M' kinks. The antikink-breather solution 
can be constructed similarly. 

For the above three types of solutions, from (2.18b) and u^ from (2.31) can be given 
explicitly in terms of the tau functions f,g and their complex conjugate as 

= In 1^, (2.37) 

{g*gf - {f*ff 

Note that (2.37) provides real solutions of equation (2.11). 

3. Properties of solutions 

In this section, we describe the properties of real solutions constructed in section 2. We 
address both the kink and breather solutions. 

3.1. 1-soliton solutions 

The tau-functions for the 1-soliton solutions are given by (2.32)and (2.33) with A^ = 1: 

f = l + ie^'+'^\ (3.1a) 

^ = l + ie«i-'^i, (3.16) 

with 

ii=Piy + - + U. e''^ = J\±^. (3.1c) 

Pi V ^^1 

The real parameters pi and .^lo are related to the amplitude and phase of the soliton, 
respectively and ^i is the phase variable characterizing the solution. The parametric 
representation of the solution (2.35) can be written in the form 

u = 2 tan^^ ( i7l + pf sinh {i J + tt, (3.2a) 

11 



X = y — T + 2 tan {pi tanh ,^1) + 2 tan pi + i/q. (3.26) 

Note that if u solves equation (1.2), then so do the functions ±u + 27rn {n : integer). For 

investigating solutions of traveling-wave type like 1-soliton solutions, it is convenient to 

parameterize solutions in terms of single variable .^1. To this end, we introduce a new 

variable X by 

X = X + Cit + Xq = 1-2 tan~^(pi tanh ^1) -|- i/q, (3.3a) 

Pi 

where 

ci = ^ + 1, (3.36) 

pi 

and xq = ^lo/pi — 2 tan~^pi. Here, we used (2.3) and (3.2b). Observing the soliton in 

the original {x,t) coordinate system, it travels to the left at the constant velocity ci. The 

soliton takes the form of a kink or an antikink depending on the sign of pi. To see this, 

we compute ux by using (3.2a) and (3.3a) to obtain 

2pi cosh ^1 
'^x = n o -. :^- (3-4) 



Vl + Pfcosh^^i + Y 



^ 



This expression implies that if pi > 0, then the solution u becomes a monotonically 

increasing function of X and has the boundary values u{—oo) = 0,-u(-|-oo) = 27r. If 

Pi < 0, on the other hand, it represents an antkink solution. Figure 1 shows a typical 

profile of the kink solution as a function of X together with the corresponding profile of 

V = Ux- To study the propagation characteristic of the soliton, we derive the dependence 

of the soliton velocity on the amplitude. To this end, let A{> 0) be the amplitude of v. 

It follows from (3.4) that 

A^WlTI^VcL, (3.5) 

2pl + l ci + 1' ^ ^ 

where, in passing to the last line, we used the relation (3.3b). Notice from (3.5) and 

ci > 1 by (3.3b) that < A < I. Solving (3.5) for ci gives 



ci = — [-A^ + 2 + 2 Vr^A2] . (3.6) 

We see from (3.6) that the velocity of the soliton is a monotonically decreasing function 
of the amplitude. In other words, the small soliton travels faster than the large soliton. 

12 



Note, however, that if one transforms to the laboratory coordinate system {Z, T) defined 
by the relations Z = x + t^T = x — t, then the velocity Ci of the soliton uz turns out to be 
a monotonically increasing function of the amplitude A. Indeed, expression corresponding 
to (3.6) becomes Ci = 1 — 2/(A^ + 1)) ^ > 1- This feature is the same as that of the sG 
soliton solution expressed in terms of the laboratory coordinate. The situation is different 
for soliton solutions of equation (1.1) with v = —1, as detailed in I. In this case, the 
velocity would become a monotonically decreasing function of the amplitude for certain 
range of the amplitude parameter. 

In view of the importance of the 1-soliton solution as an elementary solution, we 
provide an alternative derivation of the solution in appendix. The derivation is simpler 
compared with that presented in section 2 and is used frequently in reducing PDEs to 
tractable ordinary differential equations (ODEs). 

Figure 1 

3.2. 2-soliton solutions 

The tau-functions for the 2-soliton solutions read from (2.32) and (2.33) with A^ = 2 in 
the form 

/ = 1 + i (e^i+'^i + e«^+'^^) - <5e«l+«2+^l+^^ (3.7a) 

^ = 1 + i (e«i-'^i + e«2-^2) _ ^g6+6-di-d2^ (3 75) 

with 



&=ft1/+ - + «.«. e^' = jii^ = 1.2). «=j^i^. (3.7e) 

The parametric solution (2.35) with (3.7) represents three types of solutions, depending on 
values of the parameters pj and C,oj {j = 1,2), i.e., kink-kink, kink-antikink and breather 
solutions. 

3.2.1. Kink-kink solution 

If we specify pi and p2 be positive and ,^oi and ,^02 be real, then the kink-kink solution is 
obtained. The solution represents the so-called 47r kink. In figure 2a-c, we depict a typical 

13 



profile of v{= u^) instead of u for three different times. It represents tlie interaction of 
two solitons with the amphtudes Ai = 0.38 and A2 = 0.75. As evidenced from figure 2, 
a smaller soliton overtakes, interacts and emerges ahead of a larger soliton. This reflects 
the fact that the velocity of each soliton is a monotonically decreasing function of its 
amplitude (see (3.6)). The general formula for the phase shift arising from the interaction 
of A^ solitons will be given by (3.18) below. In particular, for A^ = 2, it reads 



Ai = \n[^ — ^] +4tan-ip2, (3.8a) 

Pi \Pl+P2, 



2 

I / r)i — rin \ 

A2 = 

P2 



Ll^fPl — Pl] -4tan-Vi- (3.86) 

?2 \Pl+P2j 

It can be verified from (3.8) that Ai > and A2 < for < pi < p2- In the present 
example, formula (3.8) yields Ai = 10.3 and A2 = —4.2. If one observes the interaction 
process in the laboratory coordinate system introduced in section 3.1, then one can see 
that the larger (smaller) soliton suffers a positive (negative) phase shift after the inter- 
action, which is in accordance with the property of the sG 2-soliton solution written in 
terms of the laboratory coordinate. 

Figure 2 a-c 

3.2.2. Breather solution 

The breather solution can be constructed following the parameterization given by (2.36). 
For M = 1, let 

Pi = a + 16, P2 = a — ib = pi, (a > 0, b > 0), (3.9a) 

^10 = A + i/i, ^20 = A - i/i = Cio- (3.96) 

Then, / and g from (2.32) and (2.33) become 

/ = 1 + i(e«i+'^i + e^i-'^i) + ( - J e«i+«i+"'i~'^i, (3.10a) 

14 



(3.106) 



where 



e = a{y + 



1 



a2 + 62 



r +X, 



X = b{y- 



a^ + b^ 



r +/i, 



,rfi 



'l-a2-62 + 2ia 



ae 



i/3 



(3.10c) 
(3.10rf) 

(S.lOe) 
(3.10/) 



a2 + (1 - 6)2 

The tau functions / and g can be written in terms of the new variables defined by (3.10) 
as 2 

/ = l + |-asin(x + /?) + - sm{x-f3)\e'+ (-j e^' cos 2f3 



+i 



a cos(x + /?) + - cos(x - /3) U'' + ( - 1 e^" sin 2/3 



(3.11a) 



g = l+ {a sin(x + /9) sin(x - l3)\ e^ + ( - ] e^" cos 2/3 



a 



+i 



,261 



a cos(x + /3) + - cos(x - /?) i e^ - ( - ] e^'^ sin 2/3 
a I \ a 



(3.116) 



The phase variable 6 characterizes the envelope of the breather whereas the phase variable 
X governs the internal oscillation. The parametric solution can be written by (2.35). In 
figure 3a-c, a typical profile of u is depicted for three different times. We see that the 
breather propagates to the left while changing its profile. The propagation characteristic 
of the breather is similar to that presented in I. 



Figure 3 a-c 



3.3. N-soliton solutions 



The solutions including an arbitrary number of solitons can be constracted from the 
parametric representation (2.35) with tau functions (2.32) and (2.33). There exist a 
variety of solutions which are composed of any combination of kink, antikink and breather 
solutions. Here, we address the A^-kink solutions and M breather solutions. For the former 
solutions, we investigate the asymptotic behavior of solutions for large time and derive the 
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formulas for the phase shift while for the latter ones, we provide a recipe for constructing 
M breather solution from the A^-soliton solution. As an example, we present a solution 
describing the interaction between a kink and a breather. 

3.3.1. N-kink solution 

Let the velocity of the jth kink be Cj = (l/pj) + 1 {pj > 0) and order the magnitude of 
the velocity of each kink as Ci > C2 > ... > Cat. We observe the interaction of A^ kinks 
in a moving frame with a constant velocity c„. We take the limit t — )■ —00 with the 
phase variable ^„ being fixed. We then find that / and g have the following leading-order 
asymptotics 

/ ~ (5„ exp 



jg?n+dn+5i ' 



g ~ Sn exp 



Y. (fe+rf, + f.) (1 

=n+l J 

N 1 



.jr'=n+l 



(-) 



where 



,t1, \Pn+pJ 



Sn 



j=n+l 



n 



Pj - Pk 
Pj + Pk 



(3.12a) 
(3.126) 

(3.12c) 
(3.12rf) 



n+l<j<k<N 

If we substitute (3.12) into (2.35), we obtain the asymptotic form of u and x: 



■u ~ 2 tan '^ 



T^sinh(e„ + (5i-)) 



+ 7r. 



(3.13a) 



N 



X ^ y — T + 2 tan ^ [p„ tanh (,^„ + 5^ ')] +4 N^ tan ^ pj + 2 tan ^ Pn + Vo- (3.136) 

jr'=n+l 

As t — )■ +00, the expressions corresponding to (3.13) are given by 



■u ~ 2 tan -^ 



l+p2sinh(a + 5W) 



TT, 



(3.14a) 



71-1 



X ~ y - r + 2 tan"^ [p^tanh (^„ + 5^+^)] + 4Etan"Vj + 2 tan^Vn + Vo- (3.146) 
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with 

n-l , N 2 

Let Xc be the center position of the nth kink in the (x, t) coordinate system. It simply 
stems from the relation .^n + <5n = by invoking (3.13a) and (3.14a). Thus, as t — ?■ — oo 

1 _ ^ _ 

where a;„o = ino/Vn — 2 tan^^p„. As t — )> +oo, on the other hand, the corresponding 
expression turns out to be 

-. n—l 

Xc + Cut + Xno ~ 5^'^-' + 4 ^ tan^^ pj + i/o- (3.16) 

Pn i=i 

If we take into account the fact that all kinks propagate to the left, we can define the 

phase shift of the nth kink as 

A„ = Xc(i ^ -oo) -Xc(t ^ +oo). (3.17) 

Using (3.12c), (3.14c), (3.15) and (3.16), we find that 

N n-l 

+4 ^ tan^Vj- -4 J^tan-^i, (n = 1, 2, ..., A^). (3.18) 

j=n+l j=\ 

The first term on the right-hand side of (3.18) coincides with the formula for the phase 
shift arising from the interaction of A^ kinks of the sG equation [6, 7, 9] whereas the 
second and third terms appear as a consequence of the coordinate transformation (2.3). 

3.3.2. M-breather solution 

The construction of the M-breather solution can be done following the similar procedure 
to that for the 1-breather solution developed in section 3.2.2. To proceed, we specify the 
parameters in (2.32) and (2.33) for the tau-functions / and g as 

P2J-1 = P2j = Oj + ibj, ttj > 0, bj > 0, (j = 1, 2, ..., M), (3.19a) 
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6,-1,0 = e2*,,o = A, + i/^„ (j = 1, 2, ..., M). (3.196) 

Then, the phase variables ^2j~i and C,2j are written as 

6,-1 = e, + ixj, (j = 1, 2, .., M), (3.20a) 

C2j=e,-ixj, (j = l,2,...,M), (3.206) 

with the real phase variables 

e, = aj{y + Cjt) + X„ (j = 1, 2, ..., M), (3.20c) 

Xj = b,{y - c,t) + /i,-, (j = 1, 2, ..., M), (3.20rf) 

^^• = ^T7;2' (J = 1,2,...,M). (3.20e) 

^3 + "i 

The parametric solution (2.35) with (3.19) and (3.20) describes multiple collisions 
of M breathers. One can perform an asymptotic analysis for the M-breather solution, 
showing that the M-breather solution splits into M single breathers as t — ?> ±oo. The 
resulting asymptotic form of the solution is, however, too complicated to write down and 
hence we omit the detail. One can refer to the similar analysis to that for the M-breather 
solution of the short pulse equation [10]. 

3.3.3 Kink-breather solution 

We take a 3-soliton solution with parameters pj and C,oj {j = 1,2,3). If one impose the 
conditions that p2 = p*, 62 = Qi as already specified for the breather solution (see section 
3.2.2) and P3{> 0), .^03 I'eal for the kink solution, then the expression of -u would represent 
a solution describing the interaction between a kink and a breather. We choose Pi,P2,C,io 
and 60 as those given by (3.9). Then, the tau functions / and g from (2.32) and (2.33) 
become 

/ = 1 + i (sie«^ + ^S + S3e«3^ + Q^ ' |e6+« 

-Susisse^'^^^^ - 5*3^e«+«3 + j j^^V '-1^6^s6l,e^^^^'^^^ (3.21a) 

s^ \a/ s^ 

g = l + i f le«^ + .te« + -e^A + ( -] ' ^e«^+« 
V^i S3 J \aj si 
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^^^rHs _ 5*3£le«i*+«3 +i(^Y ^<5i3<5*3e«^+«+«^ (3.216) 

'1S3 S3 \aj S1S3 



where 



1 — 6 + ia 1 1 + ips f a — p^ + ib 



2 



Si = e ' = W 7 = — , ^3 = W T —, Sis = -— —7 = '523- 

Vl-o-ia s*2 yi-ip3 \a + p3 + ibj 

(3.21c) 

Figure 4a-c shows a typical profile oi v = u^ for three different times. We see that the 

soliton overtakes the breather whereby it suffers a phase shift. An asymptotic analysis 

using the tau functions (3.21) yields the formula for the phase shift of the soliton, which 

we denote A. Actually, one has for p^ < a"^ + b"^ 



— In — — + 4 tan 



A = — Inp '- — -+4tan-^- — -, (3.22a) 

Ps [Ps ~ ci) + c^ 1 — a"^ — 0"^ 



and for a + b < p^ 



A, -Imfe + "); + ";_ 4 tan- 



2a 



(3.226) 



P3 (P3 ~ fl)^ + ^^ 1 — a^ — 6^ 

In the present example, formula (3.22a) gives A = 7.7. 

Figure 4 a-c 

3.3.4 Breather-breather solution 

The breather-breather (or 2-breather) solution is reduced from a 4-soliton solution follow- 
ing the procedure described in section 3.3.2. Figure 5a-c shows a typical profile of u for 
three different times. It represents a typical feature common to the interaction of solitons, 
i.e., each breather recovers its profile after collision. 

Figure 5a-c 
4. Reduction to the short pulse and sG equations 

We write the short pulse equation in the form 

utx = u- -^iu^)xx, (4.1) 
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where u = u{x,t) represents the magnitude of the electric field and z/ is a real constant. 
The short pulse equation (4.1) with u = —1 was proposed as a model nonlinear equation 
describing the propagation of ultra-short optical pulses in nonlinear media [11]. Quite 
recently, equation (4.1) with u = 1 

utx=u- -{u^)^^, (4.2) 

o 

was shown to model the evolution of ultra-short pulses in the band gap of nonlinear meta- 
materials [12]. Here, we demonstrate that the generalized sG equation (1.2) is reduced to 
an alternative version of the short pulse equation (4.2) by taking an appropriate scaling 
limit combined with a coordinate transformation. The reduced forms of equations corre- 
sponding to (2.7), (2.9), (2.10) and (2.11) are also presented. The iV-soliton solution of 
the short pulse equation can be derived from that of the generalized sG equation. The 
reduction to the sG equation is briefiy discussed. 

4-1. Reduction to the short pulse equation 

4.1.1. Scaling limit of the generalized sG equation 

The reduction to the short pulse equation (4.2) can be done by employing the procedure 
developed in I. Therefore, we outline the result. Let us first introduce new variables with 
bar according to the relations 

u = -, x = -{x + t), y = - yo = —, t = et, r = er, 

e e e e 

Pj=m^ Oo = 0o, (j = l,2,...,iV), (4.3) 

where e is a small parameter and the quantities with bar are assumed to be order 1. 
Rewriting equation (1.2) in terms of the new variables and expanding sin tu in an infinite 
series with respect to e and comparing terms of order e on both sides, we obtain equation 
(4.2) written by the new variables. 

Under the scaling (4.3), expression (2.7) is invariant and hence we put = to give 

Uy = sinh 0. (4.4) 

Equation (2.9) then reduces to 

(j)f = u. (4.5) 
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Equations (2.10) and (2.11) now become 



u 



ry 



VTT 



u 



uf- 



^ry 



sinh 



(4.6) 
(4.7) 



respectively. Equation (4.7) is known as the sinh-Gordon equation. 

4- 1-2. Scaling limit of the N-soliton solution 

To derive the scahng hmit of the A^-sohton solution, we use the expansion 



exp 



' N \ JV /, ip, 

Em- =n ^ 



iPj \ 2 



exp 






, , Pi 

,j = l / jr = l \ e / \ J = l / \ 3 = 1 •> . 

as well as the scaled variables (4.3). These are substituted into (2.32a) to obtain the 
expansion of the tau function / 



N 



/i=0,l \ j=l ^^ 

= /-ie/, + 0(e^), 



N 



.jr = l l<j<fc<A' 



where 



/ = ^ exp 

M=0,1 



N 






r 



^i = Pi^ + - + ^io, (j = l,2,...,Ar), 
Pi 



:.7ife 



Vj - Pk 



U,k = l,2,...,N-Jy^k). 



.Pi + Pk . 
Similarly, it follows from (2.32b), (2.32c), (2.32d) and (4.8) that 

g = g + iegf + 0{e^). 
g' = f + ief, + 0{e'), 



0(e' 



(4.9a) 

(4.96) 
(4.9c) 
(4.9rf) 

(4.10a) 

(4.106) 
(4.10c) 
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with 



9 



)_Jexp 


" N 

5^ ^^j^j + XI ^^i^^kljk 


M=0,1 


Li=i i<i<fe<JV J 



(4.10rf) 



By introducing (4.3), (4.9) and (4.10) into (2.29) and taking the hmit e — ?■ 0, we obtain 
the parametric solution of the short pulse equation (4.2) in terms of the tau functions / 
and g as follows: 

u = 2(\n^ , (4.11a) 

x = y-2{\nfg)r + yo. (4.116) 

If the tau functions are real, then (4.11) gives rise to real solutions. The construction 
of real solutions can be done by means of the procedure developed for the short pulse 
equation [10, 13, 14]. 

We exemplify real solutions by taking all parameters be real. It turns out that the 
solution develops singularities. We can confirm this fact by considering the simplest 
1-soliton solution. In fact, it follows from (4.9) and (4.10) that the corresponding tau 
functions are given by 

/ = l-e«S ^ = l+e«S ei=Piy+-+6o, (4.12) 

Pi 

where the bar appended to the variables is omitted for simplicity. The parametric solution 
(4.11) now becomes 

u = --— ^, (4.13a) 

Pi smh 4i 

2 1 2 

x = y — T-T Vyo- (4.136) 

Pi tanh ^1 pi 

Although Ux takes finite value as expected from the scaling limit of the relation (2.31), u 
itself diverges as |x| — )> oo. The general A^-soliton solution exhibits the singular nature. 
Note from (4.5) and (4.11a) that = 2hi{g/ f) gives the singular A^-soliton solution of 
the sinh-Gordon equation (4.7) [15]. The more detailed description of solutions will be 
reported elsewhere. 

4.-2. Reduction to the sG equation 
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If we introduce the following new scaled variables 



u = u, X = ex, y = ey, t 



T 

e 



Pj = ^, Oo = 0o, (j = l,2,..,iV), 



(4.14) 



then in the limit of e — )■ 0, we can deduce the generalized sG equation (1.2) to the sG 
equation 

Utx = sinu. (4-15) 

The scaling limit of (2.29b) now leads to the expression y = x which, combined with the 
obvious relation f = i, yields the limiting form of the tau functions (2.32) and (2.33) 



/ = /, f = f, 9 = f, 9 = /', 



where 



/ = 


= > ; exp 

^=0,1 


f- 


= )^exp 

M=0,1 



j=l l<j<k<N 

" N 



^<3<k<N 



t 



^j = PjX + — + ^jo, {j = 1,2,...,N), 
Pj 



^ijk 



Pj - Pk 



{j,k = l,2,...,N;j^k). 



(4.16a) 

(4.166) 

(4.16c) 

(4.16rf) 
(4.16e) 



.Pj + Pk . 

The parametric solution (2.29) with the tau functions (2.32) and (2.33) reduces to the 
usual form of the iV-soliton solution of the sG equation i.e., 



M(x,f) = 2i In ^. 



(4.17) 



As for the scaling limit of the amplitude, the appropriate scaling is A = eA. This 
relation and the scaling pi = epi from (4.14) are introduced into (3.5) to obtain the 
limiting relation A = 2pi, pi > 0, which is just the amplitude of the sG soliton in the 
{x,t) coordinate system. The phase shift is scaled by A„ = eA„. The limiting form of 
the phase shift is given by the first term on the right-hand side of (3.18), reproducing the 
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formula for the A^-soliton solution of the sG equation [6, 8, 9]. Thus, the sG limit can be 
performed consistently. 

5. Conservation laws 

The generalized sG equation (1.2) possesses an infinite number of conservation laws. Their 
construction can be done following the similar procedure to that developed for equation 
(1.1) with z/ = — 1 [2]. Here, we shall demonstrate it shortly. 
First, let 

a = u — i sinh"^ Uy. (5.1) 

By direct substitution, we find the relation 

ary - sin a = \{^ + ul)^ - i^> < UL^ - sin -u > . (5.2) 



dy j \ (1 + u. 



y' 



Thus, if -u is a solution of equation (2.10), then a given by (5.1) satisfies the sG equation 
(2.13a). Recall that equation (2.10) is a transformed form of an integrable equation (1.2) 
by means of the hodograph transformation (2.3). Thus, relation (5.1) gives a Backlund 
transformation between solutions u and a of the two integrable equations. This observa- 
tion allows us to obtain conservation laws of equation (1.2) quite simply. First, note that 
the sG equation (2.13a) admits local conservation laws of the form [16, 17] 

Pn,r = Qn,y, (u = 0, 1, 2, ...), (5.3) 

where P„ and Qn are polynomials of a and its y-derivatives. Rewriting this relation in 
terms of the original variables x and t by (2.4) and using equation (2.2), we can recast 
(5.3) to the form 

{rPn)t = {rPn cos u + Qn)x- (5.4) 



The quantities 



r"00 



— oo 



rPndx, (n = 0,l,2,...), (5.5) 

then become the conservation laws of equation (1.2) upon substitution of (5.1). The 
explicit calculation of conservation laws can be done straightforwardly. We present the 
first three of them. The corresponding P„ for the sG equation may be written as [16, 17] 

Po = l-cos(7, Pi = \^l P^ = \^t-^ly (5-6) 
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It follows from (5.5), (5.6) and the relations r^. = —u^Uxx/r, {u^ 


./r) 


r 


Uxx/r^ 


which stem 


from (2.1) that 

/»oo 










Jo = / (r — cos u)dx, 

J ~oo 








(5.7a) 


J _l r (ul ni\ 








(5.76) 


J -oo 


\lut 3ul 1 / 2 5 2 \ , 7<x 
[4r3 ' 2 r5 ' r^ V^^^ 2 "V ' r9 
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4 




dx. 


(5.7c) 



Note that terms including the imaginary unit i do not appear in (5.7) which take the form 
of total differential and become zero after integration with respect to x. 

The conservation laws generated by the procedure outlined above reduce to those of 
the short pulse and sG equations in the scaling limits described in section 4. In particular, 
the first three conservation laws of the short pulse equation (4.2) read 



(r - l)dx, (5.8a) 

^1 = -^ / ^dx, (5.86) 

dx. (5.8c) 





1 


^oo 


V^ 






^ _ 




/ 




d.x. 






2J 


f 

— oo 


r^ 






Ulr, 


- + 


lu\,. 
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6. Conclusion 

In this paper, we have developed a systematic procedure for solving the generalized sG 
equation (1.2). The structure of solutions was found to differ substantially from that of 
the generalized sG equation (1.1) with u = —1 which has been detailed in I. We have 
presented three types of solutions, i.e., kink, breather and kink-breather solutions and 
investigated their properties. We emphasize that equation (1.2) does not admit singular 
solutions in the sense that solutions must have finite slope as required by the relation 
(2.31). Consequently, loop solitons and other types of multi- valued solutions never exist. 
The existence of multisoliton solutions and an infinite number of conservation laws 
strongly support the complete integrability of the equation although its rigorous proof 
must be discussed in a different mathematical context. Another interesting issue to be 
resolved in a future work will be the initial value problem. In the case of equation (1.1) 
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with u = —1, the solution to the problem can be expressed by the solution of a matrix 
Riemann-Hilbert problem [2]. At present, however, whether the method employed in [2] 
works well or not for equation (1.2) is not known. As for solutions, the bilinear formalism 
used here and in I can be applied to equation (1.1) as well to obtain periodic solutions. 
Actually, some periodic solutions have been presented for the short pulse equation (4.19) 
with u = —1 [13, 14]. We expect that the generalized sG equation exhibits a variety of 
periodic solutions when compared with those of the sG equation. These problems are 
currently under study. 
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Appendix A. An alternative derivation of the 1-soliton solutions 

The 1-soliton solutions take the form of traveling wave 

u = u{X), X = x + Cit + xo. (A.l) 

Substituting this expression into equation (1.2) and integrating the resultant ODE once 
with respect to X under the boundary condition ^(—00) = (mod27r), we obtain 

2 _ (ci-cosm)^-(ci-I)^ 

Since u'j^ > 0, we must require that the right-hand side of (A. 2) is nonnegative. One 
can see that this condition becomes ci > cos^('u/2). In accordance with (3.3b), we solve 
equation (A. 2) under the condition ci > 1. To proceed, we define a new variable C, by 



X = / (ci - cos u)d^. {A.3) 

Then, equation (A. 2) reduces to 

u^ = ±v/(ci-cosm)2-(ci-1)2. (AA) 

Equation (A. 4) is integrated through the change of the variable s = tan(-u/2). After a 
few calculations, we obtain 



s = ±,r-^^ ^=- (A5) 

^' ci sinhv/HT^T^ ^ ^ 



and 

1 - 5=^ ^ ^ 2(ci - 1 

1 + s2 ci sinh^ a/ci - 1 ^ + Ci - 1 

Substituting (A. 6) into (A.3) and performing the integration with respect to ^, we find 



cosM=- ^ = 1 , ;^ ' . (A6) 



X = (ci - 1)^ + 2 tan"^ I tanh(v/ci-in + Vo, 



iA.7) 



where yo is an integration constant. It follows from (A. 5) and the boundary condition for 
u that 

M = 2tan~M./ ^^ sinhVci- U) +7r. {A.8) 

If we put Ci = (1/pi) + 1 (pi > 0) and C, = pi^i, then we can see that (A. 7) and (A.8) 
coincide with (3.3) and (3.2a), respectively. 
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Figure captions 

Figure 1. The profile of a kink u (solid line) and corresponding profile oiv = ux (broken line). 
The parameter pi is set to 0.4 and the parameter y^ is chosen such that the center 
position of Mx is at X = 0. 

Figure 2. The profile of a two-soliton solution v = u^ for three different times, a: t = 0, b: 
t = 2, c: t = 4. The parameters are chosen as pi = 0.2, p2 = 0.5, .^lo = —8, ,^20 = 0. 

Figure 3. The profile of a breather solution for three different times, a: t = 0, b: t = 5, c: 
t = 10. The parameters are chosen as pi = 0.3 + 0.5 i, P2 = Pi = 0.3 — 0.5 i, ,^10 = 

Figure 4. The profile oi v = u^ for three different times which represents the interaction 
between a soliton and a breather, a: t = 0, h: t = 15, c: t = 30. The parameters are 
chosen as pi = 0.2 + 0.4 i, P2 = pI = 0.2 - 0.4 i, p^ = 0.3, ^lo = 60 = 0, 60 = -30. 

Figure 5. The profile of a breather-breather solution u for three different times, a: t = 0, 
b: t = 15, c: t = 30. The parameters are chosen as pi = 0.1 + 0.2i,p2 = Pi = 
0.1 - 0.2i,p3 = 0.15 + 0.3 i, P4 = Pl = 0.15 - 0.3i, 60 = ^20 = "15, Cso = Co = 0- 
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Figure 2 a-c 
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